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Bianchi -I, -III, and FRW type models minimally coupled to a massive spatially homogeneous scalar field (i.e. 
a particle) are studied in the framework of semiclassical quantum gravity. In a first step we discuss the solutions 
of the corresponding equation for a Schrodinger particle propagating on a classical background. 
The back reaction of the Schrodinger particle on the classical metric is calculated by means of the Wigner function 
and by means of the expectation value of the energy-momentum-tensor of the field as a source. Both methods in 
general lead to different results. 
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1. Introduction 

Applying the well-known Dirac quantization 
procedure to general relativity coupled to a free 
massive scalar field leads to the corresponding 
Wheeler-DeWitt equation and to the diffeomor- 
phism constraints. Expanding the wave func- 
tional in powers of the gravitational constant 
yields a semiclassical approximation Q , which 
describes the quantized matter field in a classical 
curved spacetime, i.e. the quantized field prop- 
agate on a classical background defined by the 
gravitational degrees of freedom. It is then natu- 
ral to ask in which way the background metric is 
influenced by the quantized field. Unfortunately 
no general procedure is yet known for calculating 
this kind of back reaction from quantum gravity. 

The aim of this contribution is the following 
one: 

a) A discussion of the semiclassical approx- 
imation for different types of minisuperspaces 
especially for Bianchi -I, -III, and Friedman- 
Robertson- Walker type models minimally cou- 
pled to a massive particle. 



*H. W. thanks the DFG Graduiertenkolleg for financial 
support. 



b) To compare the results of two different ap- 
proaches for calculating the back reaction in de- 
tail: The back reaction is calculated first with 
help of Wigncr's function and second by using 
the expectation value of the energy-momentum- 
tensor of the particle as a source. 

The paper is organized as follows: 

In chapter 2 we sketch the main features of 
semiclassical quantum gravity. 

In chapter 3 the solutions of the corresponding 
Schrodinger equation of the quantized particle are 
given. 

In chapter 4 we calculate and compare the back 
reaction on the metric by the different methods 
mentioned above. 

2. Semiclassical Quantum Gravity 

Starting from the Wheeler-DeWitt equation 
for minisuperspaces 

-h- l > 2 G 2 tfG abcd d hab d Ki 

-Vh(R- 2k) + GH M ]y(hab,<j)) = 0, (1) 

where G is 167T times the gravitational constant, 
h a b the homogeneous three-metric, h its deter- 
minant, R the three-dimensional Ricci scalar, A 
the cosmological constant, G a b c d the metric on 
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superspace, dh ab = -^-^ and H M the Matter- 
Hamiltonian for the massive spatially homoge- 
neous field <j){t) 

H M = -fi 2 {2^hY x dl + Vhm 2 <t) 2 /2. (2) 
We expand the wave function \& in powers of G 

* = e*v{i[G- 1 S_ 1 (h ab ) + S {h abl <j>) + . . .]/ft},(3) 

where is a solution of the Hamilton- J acobi 
equation of pure gravity. Defining the function 
/(/i a b, 0) and the semiclassical time t(h ab ) by 

f(h ab , 4>) = F(h ab ) exp(iS (h ab , <j>)/h), (4) 

d t = G abcd {d hab S^i)d hcd , (5) 



where F{h ab ) is determined by 5_i Jlj] [g[ , we ob- 
tain the Schrodinger equation on curved space 



ihd t f{<t>,t) = H M {<l),t,d (t ,)f{ct>,t). 



(6) 



Since we are dealing with a time dependent 
Schrodinger equation, we use the Schrodinger 
inner-product for the wave function f(4>,t) 



</|/>= / #/(0,/)/*(<M), (7) 

J — OG 

where /* denotes the complex conjugated of /. 

3. Explicit solvable Models 

Our starting point is the Schrodinger equation 
(^) for a massive field. A redefinition of the semi- 
classical time d T = 2^/hdt yields 

ihdrf = [h 2 dl + 4> 2 V{t)\ /, V(r) = h{r) (8) 

with a time-depent potential V(t). Inserting the 
special ansatz 

/o = cxp (-itfdr ln[\y(r)\}/m - ln[|t/(r)|]/2), (9) 

yields for the unknown function y(r) a second 
order ordinary differential eq. 



d 2 y{r) 

d T 2 



-V(r)y(r). 



(10) 



Assuming we have found a solution of (JlC|) , we 
insert the unknown function rj(4>, r) defined by 
f(<t>,r) = f (4>,T)r)((t),T) in eq. (§). After a 



Fourier- Transformation 
elude 



p, r] —> r\ we con- 



d T fj — p(d T ln(y))d p fj — ihp 2 rj = 0. (11) 
Thus solutions of (||) are given by 

r+oo 

/ = /o / dpyexp{ip(j)} 

(12) 



exp{ihp 2 y 2 J dry 2 (f)}F(yp), 



with an arbitrary function F{yp). A complete 
set of eigenfunctions with respect to the inner- 
product (0) is given by 

f n =y- 1 / 2 H n [<P(m q (r))- 1 / 2 } 

cxp[-i(f) 2 {d T \n(y})/'ih+nJdfq- 2 (f)/2},(13) 

where H n are the Hermite polynomials of order 
n and q(r) is defined by 



q(r) =y (r) / dry (f) 



(14) 



Thus, the solutions of the Schrodinger eq. (^) 
can be calculated from the solutions y(r) of eq. 
( ^0| ) . The solutions of eq. (|l(]) have to be chosen 
in such a way that SF£(— id T ln(y)) is positive i.e. 
the functions /„ (|l3|) are normalizable. 



3.1. FRW Models 

A line element of the homogeneous and iso- 
tropic Friedman-Robertson- Walker models is 

ds 2 = N 2 (t)dt 2 - a 2 (t)[(l - fcr)" 1 dr 2 

+ r 2 (d9 2 + sin 2 (9) d<p 2 )}, (15) 

with fc = ±l for the closed and for the hyperbolic 
universes respectively. For vanishing k we obtain 
the well-known DeSitter model, which describes a 
flat three dimensional space. For nonvanishing k 
and cosmological constant A we get a complicated 
expression for the semiclassical time parameter r 



r = {k^a^y/Aa 2 - k 



+ Ak~ 3/2 &ict&n(x/ Ak^a 2 + l)}/4. (16) 

Instead of this we introduce another parameter 
x by x = kAa 2 and obtain from the eq. ( |l0| ) for 
y(r) the well-known hypergeometric differential 
equation 

x{x-l)d 2 x y+{bx/2-2)d x y + m 2 y/{AA) = 0.(17) 
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For vanishing A, k = —l, we get with r = a 2 /4 
Bessel's eq. 

9Ar 3 d 2 y + m 2 y = 0. (18) 

In the DeSitter case with fc = 0, r = a~ 3 /(6v / A) 
and A > eq. ( JTo| ) can be written as a special 
case of Bessel's eq. namely Euler's eq. 

9Ar 2 d 2 y + m 2 y = 0. (19) 

Thus for the DeSitter model the solutions of eqs. 

can be expressed by elementary func- 
tions, depending on the value of A/m 2 : 

A/m 2 < 4/9: y = VTr- ik °, (20) 

A/m 2 > 4/9: y = ^(fcir fc ° + ifc 2 ^ fe °), (21) 

A/m 2 = 4/9 : y = 1h(t) + ik 2 ), (22) 

with fc = v/|m 2 /(9A)-l/4|, fci,fc 2 € # and 
fcifc 2 > because of normalizability. 



3.2. Anisotropic Minisuperspaces 

Contrary to the FRW models discussed above 
each of the following two models possesses two 
gravitational variables z(t), b(t). 

3.2.1. Bianchi —I with rotational symmetry 

The line element of the Bianchi-I model with 
rotational symmetry is 

ds 2 = N 2 (t) dt 2 -z 2 (t) dr 2 - b 2 (t) [dO 2 +6 2 dcf 2 ] . (23) 

Here we obtain for the semiclassical time pa- 
rameter T 



2c 



3c 2 , z 2 ) 



In 



AAb 2 
34^ 



(24) 



&(4A6 2 - „. 

where Co is the separation constant of the solution 

S-i = c z 2 b + A6 3 /(12c ) (25) 

of the Hamilton-Jacobi eq. . Inserting r of ( pif ) in 
eq. ( |l"o| ) we obtain 

G) 2 j/ + Cl sinh- 2 (c 2 T)j/-0, (26) 

with a = 3cQm 2 / 2 /(32A) and c 2 = 3c //4, where 
I — 4:Ab 3 / (3cq)— z 2 b has to be treated as a constant. 

Eq. (|2^) is related to a hypergeometric one, 
which can be seen by a change of variables w — 
l+exp(2c 2 r). This transforms eq. ( |26| ) into 

w 2 (l - w)d 2 w y - w 2 d w y + c iy /c 2 = 0. (27) 

Defining the function 8 by 9{w) = y(w)w~ c with 
c= l/2+-\/l/4-(-ci/c 2 , 9 fulfills the hypergeomet- 
ric differential eq.. 



3.2.2. Bianchi-III with rotational symmetry 

Similar to the line element ( p3| ) above, we have 

ds 2 =Ar 2 (t)^ 2 -z 2 (i)dr 2 -fe 2 (t)[d6' 2 +sinh 2 (6')d< ? !) 2 ].(28) 

Introducing the parameter 
z = (4 C 2 z 2 + l)/(4 C 2 z 2 -l), (29) 
where C3 denotes the separation constant of 
5_i =c 3 i 2 6 + V(4c 3 ), (30) 
eq. ([n]) yields 

d x (x 2 -l)d x y + Ci (l-x) 2 y = 0, (31) 

with C4 = 2e 2 m 2 (4c 2 z 2 — 1)6 as a constant. 

This eq. ( pl| ) is known from the scattering prob- 
lem for the hydrogen molecule ion first solved 
byJaffep. 



The solutions of the Schrodinger equation (|1 
can be normalized for each of the models above. 

4. Back Reaction 

Two different approaches for calculating the 
back reaction of the quantized matter field on the 
classical metric are investigated. First, using the 
expectation value of the energy-momentum ten- 
sor as a source and second, with help of Wigner's 
function Q]. We calculate and compare the re- 
sults for different solutions of the Schrodinger eq. 
d) for the DeSitter model. 

4.1. Using the energy- momentum tensor 

One of the most familiar definitions of back re- 
action is to replace Hm by < f\Hia\f > in the 
classical Hamiltonian constraint, where f is the 
wave function (||). This is equivalent to define 



P = G- 1 d a S^ 1 + < f\d a (3(a, 0)|/ >, 



(32) 



for the gravitational variable a of the DeSitter 
model, where P denotes the classical momentum 
conjugate to a and /3(a, <fi) the phase of the wave 
function /. 

Inserting the unperturbed classical expression 
P = -G^dta 2 (N = 1) in eq. (||) we get a 
first order differential eq. for the perturbed met- 
ric a(t). Since the perturbation is exact up to 
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order G only, we set a(t) = a(t) + Ga\(t), where 
a = exp[v / A(t+to)/2] is the unperturbed DeSitter 
metric. 

Given the wave function r(a)) we can cal- 
culate ai(t) explicitly. For the " groundstate" 
fa=y~ 1 / 2 exp(— i(/) 2 d T \n(y)/(4:h)) we obtain 

A/to 2 < 4/9 : a x = -hk 2 a - 2 [k + l/(4fc )], (33) 
A/to 2 > 4/9 : ai = -^a^ 2 [a 6fco /c 2 + a- 6fe °A: 2 ], (34) 
A/to 2 = 4/9 : a x = -7Lcr 2 [fcf(ln(a)) 2 

+ fe 9 ln(a) + fe 10 ] ! (35) 

with real constants k5,ke,kT,ks,kg,k\o. In the 
massless case we have: oi = —Tik^a^ 5 , fen £ i?. 
For large universes, a — > oo, the perturbation is 
negligible ai <C a, a — > a. 

If / contains the Hermite polynomial H n we 
get 

oi = -hkla- 2 {2n + l){k + n!2"/(4fc )} (36) 

for A/to 2 < 4/9. Thus, the result @ for the 
"groundstate" is only slightly modified. All the 
perturbations © , © , © , © vanish in the 
"classical limit" % —* 0. 

To obtain the additional classical part of the 
back reaction of order G, we take the wave func- 
tion 

/ = y- 1 ' 2 exp[-^ 2 a r ln(y)/(4ft) + io0/(%) 



ic 2 



dry- 2 (r)/?i], 



(37) 



and get 

oi cla 6k °- 2 ^la- 2 -cla- 6k «- 2 -Mla 6k «- 2 p8) 

with real constants Ci, 02,03,64. Here the quan- 
tum corrections are given for large a. Multiplying 
the wave function (|37j ) by Hermite polynomials 
changes the coefficient c 2 of the quantum correc- 
tion in oi only. 

4.2. Wigner's function 

Wigner's function Fw is a generalization of a 
classical correlation function on phase space Qj. 
F\y(a,P,4>,P<p) depending on the gravitational, 
the matter variable and their momenta is defined 
by 



F w =Gmr l j dudv f*(a-Ghu/2,<f>-hv/(2m)) 

J —OO 

e- i( - PuG+p f v/m) f(a + huG/2 7 (t)+nv/{2m)).(39) 



Expanding Fw in powers of G leads in order G° to 
Fw ~S(P— G~ 1 d a S-\) from which we obtain the 
unperturbed classical expression P= — G~ 1 d a S-i 
^ . Integrating over <f>, P^ we obtain in order G 

POO 

F G = (daS^y 1 #|/| 2 5(P-G- 1 a a 5_ 1 -9 a /?),(40) 



where j3(a, <fi) is the phase of /. The peaks of Fq 
yield the relation between a and P from which we 
obtain a\(t) by integration as in the case of the 
energy-momentum-tensor . 

For the "groundstate" we calculate 

A/m 2 < 4/9 : ai = -hk 2 k Q a~ 2 , 

A/to 2 > 4/9 : a x ~ -ft a - 2 (fc|a +6fc ° -/c 2 1 a" 6fc °), 



A/to 2 = 4/9 : a x hk\ 



fcu.fcia € #-(41) 



The last two expressions are given in leading or- 
der a^> 1 only. The case A/to 2 < 4/9 yields the 
same result as (p3|). A/to 2 > 4/9 agrees in lead- 
ing order with (p4[). Only A/to 2 = 4/9 leads to a 
different back reaction compared with ( |35| ) . 

The quantum corrections induced by the wave 
function (|37j ) remain unchanged. 
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